Single Electron Transport in electrically tunable nanomagnets 
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We study a single electron transistor (SET) based upon a II- VI semiconductor quantum dot doped 
with a single Mn ion. We present evidence that this system behaves like a quantum nanomagnet 
whose total spin and magnetic anisotropy depend dramatically both on the number of carriers 
and their orbital nature. Thereby, the magnetic properties of the nanomagnet can be controlled 
electrically. Conversely, the electrical properties of this SET depend on the quantum state of the 
Mn spin, giving rise to spin-dependent charging energies and hysteresis in the Coulomb blockade 
oscillations of the linear conductance. 

PACS numbers: 73.23.Hk,85.75.-d,78.67Hc,78.55.Et 
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Nanomagnets attract interest both because of their 
intriguing behavior as relatively macroscopic quantum 
objects^ and their potential technological applications as 
magnetic bits^ and qbits^. The two fundamental prop- 
erties of a nanomagnet are the net spin of its ground 
state, S and its magnetic energy anisotropy tensor, K, 
that governs the stability of the magnetization with re- 
spect to quantum and thermal fluctuations. Although 
recent experiments show that single molecule magnets 
like Mni2'*'^ or metallic Co® nanoparticles can be probed 
in single-molecule transistor measurements, their proper- 
ties can hardly be tuned once they are fabricated. Here 
we show that a single electron transistor (SET) consist- 
ing of a II- VI quantum dot doped with a single Mn atom 
behaves like a tunable nanomagnet whose magnetization 
and anisotropy axis can be reversibly manipulated elec- 
trically. Conversely, the conductance and charging en- 
ergy of the tunable nanomagnet depend on the quantum 
state of the Mn spin and are not uniquely determined by 
the gate and the bias voltage. 

Our proposal is based on two independent progress in 
nanofabrication. On one side, the fabrication and optical 
probing of single CdTe quantum dots doped with a single 
Mn aton-J-'^. In the absence of carriers, the spin S = 5/2 
of the Mn is free. Optical excitation of electron- hole pairs 
into the dot shows that the Mn spin is exchange coupled 
to both the electron and the hole^i^iiSiiS. On the other 
side, the control of the charge state of II- VI semiconduc- 
tor quantum dots with single electron accuracy has been 
experimentally demonstrate d^^i^^ as well as in the case 
of single-Mn doped quantum dotsi^ and Mn-doped GaAs 
islandsi^. 

Hamiltonian. We consider a CdTe quantum dot (QD) 
doped with a single Mn, weakly coupled to two metal- 
lic and non-magnetic electrodes. The dot can be gated 
so that either the valence band or the conduction band 
is in resonance with the metallic reservoir and the num- 
ber of either electrons or holes is varied at will. The total 
Hamiltonian reads: H = Hqd+'Hc+'Hl+'Hr+Vl+Vb.- 
Here H-qd is the Hamiltonian for the diluted magnetic 
semiconductor (DMS) quantum dot. In analogy with the 
the standard model^^ for bulk DMS, H-qd describes con- 



fined conduction band electrons and valence holes inter- 
acting with a localized Mn spin 5 = |, denoted as M, via 
a local exchange interaction. QD carriers occupy local- 
ized spin orbitals (j)a with energy Eq which are described 
in the envelope function k ■ p approacbiSiiSiiil. In the 
case of valence band holes the 6 band Kohn-Luttinger 
Hamiltonian, including spin orbit interaction, is used as 
a starting point to build the quantum dot statesii. The 
second quantization Hamiltonian of the isolated dot de- 
scribes the states of reads: 

'^QD = 2_^ i^aSa,a' + Ja,a' M ■ Sa^a' j fafa' (1) 



Here /^ creates a band carrier in the a single particle 
state of the quantum dot, which can be either a va- 
lence band or a conduction band state. The first term 
in the Hamiltonian describes non interacting carriers in 
the dot and the second term describes the exchange cou- 
pling of the carriers and the Mn. We neglect intcrband 
exchange so that Ja,a' = Je iJa,a' = Jh) if both a and 
a' belong to the conduction band states (valence band 
states). In contrast, we include exchange processes by 
which a carrier is scattered between two different levels 
of the dot that belong to the same band. The matrix 
elements of both valence and conduction band spin den- 
sity, evaluated at the location of the Mn atom, are given 
by Sa,a'- They depend strongly on the orbital nature of 
the single particle level in question. In the case of con- 
duction band we neglect spin orbit interactions so that 
Sa.a' is rotationally invarianti^. In contrast, strong spin 
orbit interaction of the valence band makes the Mn-hole 
interaction strongly anisotropioiSiiJiiil and it varies be- 
tween different dot levels. Following previous workiSii^iii 
confinement is described by a hard wall cubic potential 
with Lz < Lx — Ly. Although real dots are not cubic, 
this simple model^*''^^ provides an excellent description 
of the Hamiltonian of the Mn spin coupled to the carriers, 
which is able to account for the non-trivial single-exciton 
PL spectra both for neutral^i^i^ and charged^'^ single-Mn 
doped CdTe QD. 

Coulomb repulsion between carriers is described within 
the constant interaction model^*: Tic = wjiQ + ^g^g + 
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Ci 



MR^2 



)2, where C = Cl + Cj^ 



Cg is the total 



capacitance to the external circuit, C^ and Cn are the ca- 
pacitances of the left and right junctions (eVe = Ml — M-R 
is the bias voltage) and Cg is the capacitance to the 
gate (with voltage Vg). Q is the extra charge in the 
dot. We do not consider dots with orbital degeneracy 
for which Coulomb correlations, neglected in this pa- 
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J2ap^p^P<yK<^ describe the metalhc dee- 



per, are relevantA2,i^. Finally, Hl 

and Hr 

trodes and Vl = J2a.k,a y<y,k,aflaka- + h.c and V^ = 

'^a p a ^(j.k.afibpa + h.c are the standard spin-conserving 

tunneling Hamiltonian that couple the metallic reservoirs 

and the dot. 

We first discuss the properties of the eigenstates \N) of 
Hqd for isolated dots (Vl = Vr = 0) with a given num- 
ber of carriers, interacting with the Mn atom. We show 
resuhs for two dots of CdTe with L^ = 60A, L^, = 80l 
and different Ly = 80 A (dot A) and Ly = 75 A (dot B), 
both doped with 1 Mn atom. The neutral dot has 6 
degenerate states, corresponding to the (25* -I- 1) equiv- 
alent spin orientations of the S — 5/2 Mn spin. This 
degeneracy is lifted in the presence of either electrons or 
holes. We focus on dots with a odd number of carriers 
(open shells) for which the interactions are stronger-i^ii^ 
and study how the magnetic anisotropy varies with the 
number of carriers. The spectra of dots with 1 electron, 
1 hole and 3 holes are shown in figs, la, lb and Ic re- 
spectively for dot A (and also QD B for the case of 1 
hole). The effect of intra- level exchange is magnified in 
the inset of fig. la. In figure lb and Ic we only show 
the low energy manifold for dot A with 1 and 3 holes. 
The 12 states of the low energy manifold for Q = ±1 
and Q = -1-3 are formed mainly by the two lowest energy 
electronic configurations of the dot with a single unpaired 
fermion coupled to the 6 Mn states. The low energy sec- 
tor of Hqd can be described by an intra-level effective 
Hamiltoniania.: Heff = fr^W + jyvVMy + fr'^W 
where t° are the Pauli matrices operating on the isospin 
space defined by the lowest energy single particle doublet. 

Both the absolute and the relative values of jx , jy and 
jz depend mostly on the spin properties of the external 
shell of the quantum dot: either conduction band level 
(Q = -1), heavy hole (Q = +\) or light hole {Q = -H3). 
Thereby, the jx.y,z can be controlled reversibly by means 
of the gate voltage in the same device. The effective 
Hamiltonian of the Mn coupled to the " master fermion" 
in dot A goes from ferromagnetic Heisenberg [jx = jy = 
jz < 0) when (Q = — 1) to antiferromagnetic Ising jx = 
Jv = 0,jz > when (Q = +1) to XXZ {jx == jy > jz)- In 
dot B similar results are obtained, with jx ^ jy, which 
provides a spin-flip term in the Q = -1-1 case, absent in 
dot A. 

Correspondingly, the spin properties of the ground 
state also change as a function Q. In the case of Q = — 1 
the Mn spin 5/2 and the conduction band electron couple 
ferromagnetically to yield a septuplet with St = 3. In 
the case of Q = -1-1 the ground state doublet corresponds 
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FIG. 1; (Color online). Low energy spectra of TIqd for dot A 
(circles) with Q = -1 (a),Q = -1-1 (b) and Q = +3 (c). Notice 
the different vertical scale in (a) and (b) In panel (b) we also 
show the spectrum of QD B with Q — +1 (squares). The 
Kohn-Luttinger parameter for CdTe are 71 — 4.14, 72 — 1.09 
and 73 = 1.62 and the spin orbit interaction is A = 950 meV. 
We take Je = -15 eV A^ and A = 60eV l^ 



to the Mn spin maximally polarized against the heavy 
hole spin, \Mz = -5/2, t) and \Mz ^ +5/2, i), both for 
dots A and B. In dot A the rest of the low-energy sector 
is formed by 5 doublets eigenstates of both Mz and r^ , 
whereas in dot B the small spin-flip interaction mixes 
\Mz = +1/2, i) and \Mz = -1/2, t). In the case with 
Q = +3, the ground state is not degenerate and the Mn 
spin is polarized in the xy plane, minimizing Mz- These 
differences reflect the spin properties of conduction band 
electron, heavy hole and light hole respectively. 

We now address how these remarkably different mag- 
netic properties occurring in the same dot are reflected 
in the electrical behavior of the SET. In analogy with 
previous wor k^^i^^'^'^1^^1^^ , we derive a quantum master 
equation for the dissipative dynamics of the reduced den- 
sity matrix pNAiit) written in the basis of many-body 
states |iV). Importantly, this quantum master equation 
includes the combined dynamics of both populations and 
coherences. The latter are important because of the in- 
trinsic many-body degeneracies of the QD spectra shown 
in fig. 1 . Assuming that the quantum dot is weakly cou- 
pled to the electronic reservoirs (sequential tunneling), 
the dissipative dynamics of the density matrix is gov- 
erned by a Markovian kernel, p{t) — Ap{t), where p can 
be casted as a vector containing both populations and co- 
herence terms. The matrix A contains information about 
dissipative dynamics of p which is governed by the rates: 
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is the Fermi function of 



Here /+ EE /t, /- EE /< 

reservoir r and n^, —1-n. 

states M with charge Q are connected with states N with 

charge Q ± 1 . . The coupling to the leads is parametrized 



The notation F^ jy,j implies that 



by Tl,r 
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VlM Nl,r, where Nl.,r is the DOS of 



the metallic reservoir. Once we obtain the steady state 
density matrix p (namely, Ap = 0), we can compute the 
average charge, magnetization and current. To lowest 
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order in Tl,r, the most general expression for the current 
can be written as / = ^^~^^ with: 



h/R = eVL/R J2 EE'5iv,W'{4/i^(^^ - EN'){N'\f^\M){M\fi\N) - nl^^iE^ - EM){N'\fi\M){M\U\N)} (3) 
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Notice that Eq. ^ includes both diagonal and non-diagonal terms in the density matrix. The latter are important 
when two degenerate states with Q are coupled to the same state of Q ± 1 via a single tunneling event. 



The steady state of a standard SET is uniquely charac- 
terized by external voltagesi^. For instance, a new charge 
is accommodated in the dot at precise values of the gate 
voltage, when the electrochemical potential of the dot 
p{N) (the energy required for adding the A^th electron 
to the dot) falls within the bias window p.^ > p{N) > pR. 
When this condition is met, the number of electrons can 
vary between iV — 1 and N resulting in a single-electron 
tunneling current. Importantly, in our case the charge 
and the conductance of the SET depend also on the quan- 
tum state of the Mn spin. 

In figure 2 we show linear conductance Go(Vg), av- 
erage charge and diagonal terms of the steady state p, 
as the gate produces the transition between charge zero 
and chage ±1 for both electrons (left panels) and holes 
(right panels). The initial Vg for the charging simula- 
tions is chosen so that only the Q = states are oc- 
cupied. This initial condition is described by a ther- 
mal p with 6 equally populated Mn spin states, Mz = 
±5/2, ±3/2, ±1/2. We ramp the gate and solve the mas- 
ter equation to obtain the steady state p, which is used 
as initial condition for the next run with higher Vg- In 
the case of electrons we obtain standard results: a single 
peak in the Go(Vg) curve occurs as the gate is ramped 
to change the charge of the dot by one unit. In fig. 2c we 
show the evolution of the steady state populations: The 
Q Q — spin states are relaxed altogether in favor of the 
7 states of the Q — 1 and St — 3 states. 

The results for holes in QD A are remarkably differ- 
ent: in the process of injection of 1 hole the Go(Vg) 
curve shows 3 peaks instead of 1. This results from the 
lack of Mn spin relaxation [M^ is a conserved quantity 
for the entire Hamiltonian including tunneling), which 
makes the steady p different from to the thermal p. As 
the gate brings down into resonance the 2 Q = +1 ground 
states (+5/2,1) and (— 5/2,|) with the 6 ground states 
with Q = 0, population tranfer only affects states with 
\Mz\ ^ 5/2 in both charge sectors (first peak in fig. 2d 
at Vg — 4.3mcV). Further increase of the gate brings the 
energy of the the Q = +1 ground state doublet below the 
2 = states with \AIz\ ^ 5/2, which are not depleted 
because M^ is conserved (fig. 2f). The population trans- 
fer only occurs when the second and third doublet of the 
g = -1-1 spectrum, with \Mz\ = 3/2 and \Mz\ = 1/2 are 



brought down in resonance with the Q = states. This 
accounts for the other two peaks in the Go{Vg) curve as 
the charge of the dot approaches +1. Hence, the charg- 
ing energy for holes depends on the absolute value of the 
spin of the Mn. The discharge simulation is done anal- 
ogously. If the initial Vg is such that there is one hole 
in the QD, the Ising interaction removes the degeneracy 
among states with different \Mz\. Only the doubly de- 
generate ground state of the Q = +1 sector is occupied in 
thermal equilibrium. As the gate is ramped to discharge 
the dot, a single peak in the conductance is obtained, 
corresponding to the resonance condition with between 
the Q = and Q = +1 states with \Mz\ = 5/2 . 

The difference between electrons and holes arises from 
the different value of a crucial time scale in the magnetic 
sigle electron transistor: the Mn spin relaxation time, Ti. 
In the case of QD A, Ti is infinite for holes (Ising cou- 
pling) which makes the steady state different from the 
thermal state. In the case of electrons (left panels) the 
transverse spin interactions make Ti comparable to the 
charge relaxation time (F^ ^) so that steady and thermal 
p are identical. In real dots doped with one hole, Ti may 
be long but not infinite. Two independent mechanisms. 
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FIG. 2: (Color online). Go(Vg) (upper panels), charge (mid- 
dle panels) and diagonal terms of the p (lower panels) for QD 
^ as a function of Vg around the Q = —1 ^^0 transition 
(left) and Q — ^-^ +1 (right). Lower panels: solid (dashed) 
lines correspond to Q = (|Q| = 1) states. Results obtained 
with Fl = Fh = O.OlmeV and ksT =0.05 meV. 
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FIG. 3: (Color online). Current and differential conductance 
as a function of bias for QD A (a,b,c) at various charge states 
and QD B (d) (see text). 



missing in the simulations shown in the right panels of 
fig. 2., yield a finite Ti. First, the Mn Ti due to super- 
exchange with other Mn spins which scales exponentially 
with the Mn density^S,. For bulk Cdo.ggsMno.oosTe we 
have Ti — I0~^s, which is a lower limit estimate for Ti 
of the QD with a single Mn. The second mechanism is the 
small^"^ transverse spin interaction, which is proportional 
to the light-hole heavy-hole mixing. We have simulated 
QD B, for which spin- flip interaction between the hole 
and the Mn is small but non-zero resulting in a finite Ti . 
If we integrate the master equation for F"^ << t << Ti 
the Go{Vg) curve displays 2 peaks and hysteretic be- 
haviour. In contrast, if we integrate the master equation 
for F~^ << Ti << t, the system reaches the equilib- 
rium state for each value of Vq so that the Go(Vg) curve 
has a single peak. Therefore, we claim that effects re- 
lated to incomplete spin relaxation of the dot will be 
observed subject to two conditions: the finite bandwith 
of the measurements should be larger than l/Ti (see for 
instance ref. (@)) and the pace at which Vq is ramped 
should be faster than Ti. 



The finite bias conductance of the device also depends 
strongly on the charge state of the dot. In figure 3 we 
show I{Vb) and ^^7- curves for dot A corresponding to 
bias-assisted single electron fluctuations between Q — 
and Q — —1 (fig 3a), Q — and Q = +1 (fig. 3b), and 
Q ~ 2 and Q = +3 (fig 3c). Fig. 3d is the analogous 
of 3b for dot B. Current flows whenever the addition of 
a fermion is permitted by energy conservation and spin 
selection rules. The former provides a link between the 
-0- curve and the energy spectra of dots shown in fig. 
1 (since the spectra of dots with Q — and Q — 2 are 
flat). Interestingly, the -^M- for electron tunneling (3- 
a) shows a zero-magnetic field splitting related to recent 
experimental observations^^. In turn, --^- has 6(7) peaks 
in fig. 3b (3d), are similar to the experimental single 
exciton PL spectrap^i^i^ii^. 

In summary, we have shown some of the equilibrium 
and non-equilibrium properties of a semiconductor quan- 
tum dot doped with a single Mn atom and wired as a 
single electron transistor. The different orbital nature 
of the conduction band electrons, heavy holes and light 
holes determines both the total spin and the magnetic 
anisotropy of the dot. In the case of holes, for which Mn 
spin flip processes are heavily inhibited, we predict differ- 
ent results for the Go(Vg) curves depending on whether 
the system is relaxed to equilibrium or not. In the case 
of the latter, we predict hysteretic Coulomb blockade os- 
cillations related to the the quantum state of the Mn 
spin. Because most of the transport properties discussed 
above are inherent to nanomagnets with long spin relax- 
ation time, our findings might be very general and have 
implications in recent experiments^i^ii^. 
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